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INTRODUCTION 


Detection 


Evernden  (1969)  examined  records  from  all  underground  explosions  before 
1966  and  detected  long-period  shear  waves  from  only  two  explosions  at  one 
station  each.  This  led  him  to  suggest  that  for  fixed  m^  the  S wave  amplitudes 
for  earthquakes  would  be  about  a factor  of  ten  greater  than  for  explosions. 

Nuttli  (1969)  reported  131  measurements  of  long-period  S wave  arrival 
time,  first  peak  amplitude,  sign,  and  period  for  HALFBEAK  and  GREELEY  at  43 
WWSSN  and  LRSM  stations.  Measurements  were  made  on  all  three  components  when 
possible.  He  was  able  to  establish  a good  distance-amplitude  relation  for 
the  first  peak  of  SH  motion  from  GREELEY  measurements  alone,  and  established 
that  there  was  a node  in  the  SH  radiation  pattern  of  HALFBEAK  and  GREELEY. 

Von  Seggern  (1972)  reviewed  the  Seismic  Data  Laboratory  (SDL)  shot 
reports  and  Geotech  earthquake  bulletins  for  the  decade  1960-1970  for  the 
LRSM  network  and  for  the  VELA  observatories.  From  these  he  was  able  to  con- 
struct histograms  of  the  amplitude  ratios  for  earthquakes  and  explosions  of 
short-period  S to  short-period  P,  long-period  S to  short-period  F,  long- 
period  S to  long-period  P,  long-period  S to  Rayleigh  wave,  and  Love  wave  to 
Rayleigh  wave.  He  found  only  8,  11,  and  9 observations  of  short  period  S, 
long-period  S and  long-period  P respectively  for  explosions  out  of  the 
explosions  studied  in  the  SDL  shot  reports,  a total  of  approximately  40 
explosions . 

The  difference  of  the  means  of  von  Seggern's  long-period  S to  short- 
period  P ratios  were  in  general  agreement  with  Evernden's,  although  there 
was  some  overlap  of  individual  ratios.  Von  Seggern  also  showed  good  dis- 
crimination using  the  ratio  of  long-period  S to  LR.  All  of  von  Seggern's 
ratios  were  uncorrected  for  distance  or  magnitude. 

Von  Seggern  and  Lambert  (1972)  reported  several  S wave  observations 
from  the  explosion  MTLROW,  and  Johns  (1971)  detected  and  measured  long-period 
S waves  from  CANNIKIN  at  6 out  of  a possible  8 LRSM  portable  stations. 
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Strauss  (1973)  determined  the  long-period  S thr.  u . 
quakes  in  Kamchatka  and  the  Kuril  Isla  d h k ^ ^ ^ f°r 

- components  a„i  :ming  the  verticai' 

amplitude.  He  plotted  q , race  Wlth  the  largest 

— - — r:  ::vir:r:i;  r four  ~d — 
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with  other  observations  With  th • 6 °W*  dlSagrees  seriously 

the  — ^ ^r  explosions  ^ ^ ^ ^ 

of  A/T  for  earthquakes.  b W the  l0West  value 


Distance-Amplitude  Relatione 

Gutenberg  (1954b,  1945c),  and  Gutenberg  and  Richter  (1956)  / 

Richter  (1958)  from  which  Figure  1 is  copied^  d n . ’ ^ 

curves  for  S as  well  as  P The  d ’ ^ dlstance-amPlitude 

, he  data  were  normalized  to  M from  LR  as  d • 

cussed  by  Gutenberg  (1945a)  Th*  da,  s dls_ 

s v.  aaa j.  lhe  data  used  were  ]9?T-iQaQ  _ .. 

Gutenberg,  1945c')  Tho  • „ station  bulletins 

g,  i^Abc),  The  instruments  in  use  at  that  nm* 

not  clear  what  th*  perIod  KOuld  hav~;slted  and  u is 
personal  communication).  Thus  it  u (Stauder, 

' tnus  it  might  be  surprising  if  Ci.*-0nkc  . 
were  approximately  valid  for  S wa  & " Utenber8  s curves 

mentation.  Howeulr,  we shall III"  “ """"  l0"“  ‘”'~ 

fot  this  application.  Evernden  (1969)  p “T  aPPe3r  “ ^ ^ 

for  long-period  s whIch  , ____  a Stance-amplitude  formula 


for  long-period  s which  a “ Stance-amplitude  formu 

K period  S which  decays  much  more  slowly  than  c „ 

Richter,  i.e.  A/T  proportional  to  a"7.  ButtJ  " GUtCnber8  “d 

ttli  (1969)  plotted  amplitude 
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surface  focus,  and  are  in  striking  disagreement  with  a a’’7  ia 
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in  this  study  we  report  a distance-amplitude  relation  for  S waves  from 
earthquakes  which  is  in  excellent  agreement  with  the  shape  of  Nutrli's 
GREELEY  data  and  which  disagrees  with  Gutenberg  and  Richter's  only  around 
A = 20°,  where  the  shape  of  our  curve  shows  the  sharp  maximum  found  by 
Nut L 1 i . Gutenberg  and  Richter's  curve  is  about  0.5  magnitude  units  smaller 
than  Nuttli's  and  our  data  around  20°. 
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Source  of  S Waves  from  Explosions 

While  SV  radiation  from  explosions  could  be  created  by  P to  S conversion 
at  the  tree  surface  and  other  horizontal  interfaces,  early  workers  in  seismic 
discrimination  had  hopes  that  SH  waves  would  not  be  present  in  explosion 
signals  and  would  prove  to  be  a powerful  diagnostic. 

In  1962  Wright  and  Carpenter  suggested  that  common  field  experience 
with  small  explosives  indicated  that  SH  waves  were  generally  observed.  To 
help  explain  this  they  exploded  small  charges  in  soft  clay  and  made  plaster 
casts  of  the  resulting  holes.  Figure  4 shows  cross-sections  of  the  plaster 
cast  as  drav/n  by  Wright  and  Carpenter.  We  see  in  the  left-hand  polar  dia- 
gram a substantial  asymmetry  about  the  shot  point.  One  would  expect  that 
this  would  be  accompanied  by  substantial  SH  motion.  Wright  and  Carpenter 
tentatively  explained  the  asymmetry  as  being  a result  of  Taylor  instability 
which  occurs  when  a light  material  is  accelerating  a heavy  one.  This 
experiment  and  accompanying  explanation  might  lead  one  to  expect  more  SH 
motion  from  explosions  in  weak  material  than  from  explosions  in  strong 
material. 

Kisslinger,  Mateker,  and  McEvilly  (1961)  made  careful  measurements  of 
the  SH  radiation  from  small  explosions  in  soil.  In  Figure  5 we  see  one  of 
their  observed  radiation  patterns,  and  a perfect  fi  (after  polarity  rever- 
sal) of  SH  waveforms  at  two  stations  144°  apart.  Kisslinger  et  al  reached 
the  general  conclusion  that,  surprisingly,  "the  radiation  pattern  seems  to 
be  tied  into  the  geography  of  the  test  site". 

The  most  difficult  set  of  their  experiments  to  understand  was  the 
finding  of  nearly  identical  SH  radiation  from  a series  of  three  shots  in 
which  the  first  0.5  lb  charge  was  tamped  in  soil  at  a shallow  depth  which 
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lcd  to  cratering.  The  hole  was  then  dug  out  into  a smooth  hemisphere  and 
tilled  with  water  in  which  another  0.5  lb  charge  was  detonated;  the  hole  was 
filled  again,  and  an  0.5  lb  charge  again  detonated. 

This  experiment  would  seem  to  rule  out  cracking  as  a possible  generation 
mechanism,  leaving  only  Taylor  instability.  The  most  obvious  reason  we  can 
imagine  for  having  a preferred  orientation  for  the  SH  radiation  resulting 
from  Taylor  instability  in  these  experiments  is  a consistent  asymmetry  in  the 
orientation  of  the  actual  explosive.  Such  an  asymmetry  might  preferentially 
digger  one  particular  mode  of  Taylor  instability.  However,  Carl  Kisslinger 
(personal  communication)  says  that  the  dynamite  sticks  were  exploded  one  on 
top  of  the  other,  not  side  by  side;  and  that  the  firing  cap  was  placed  in  the 
end  of  one  of  the  sticks  of  dynamite.  A remaining  possibility  is  that  the 
soil  is  anisotropic,  leading  to  a preferred  form  of  Taylor  instability  when 
the  water  pressure  wave  hits  the  soil  surface. 

White  and  Sengbush  (1963)  also  detected  SH  motion  from  small  explosions. 
The  data  in  Kisslinger  et  al.  (1961)  were  all  recorded  at  the  surface,  and 
as  a result  were  predominantly  of  the  Love  wave  type.  They  found  "Love  to 
Rayleigh"  ratios  on  the  order  of  1:4.  Geyer  and  Martner  (1969)  on  the  other 
hand  used  down-hole  recorders  and  thus  were  able  to  measure  SH  and  SV  directly. 
They  obtained  SH  to  SV  ratios  as  great  as  5 or  10  to  1 in  some  cases  (the 
apparent  relative  amplitudes  in  their  Figures  are  accurate;  S.  T.  Martner, 
personal  communication)  and  they  conclude: 

"1.  SH  waves  are  commonplace.  In  fact,  they  have  been  identified 
in  almost  every  area  where  we  have  used  three— component 
seismometers . 

2.  The  most  striking  development  of  the  SH  waves  is 
associated  with  horizontal,  high-contrast  interfaces. 

(In  fact,  it  might  be  suggested  that  this  could  be  a 
way  to  shoot  nuclear  devices  to  generate  large  amounts 

of  SH  waves  so  that  they  could  be  disguised  as  earthquakes.) 

3.  The  amplitude  of  SH  waves  increases  approximately 
linearly  as  a function  of  charge  weight. 
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4.  Azimuthal  data,  though  limited  in  these  experiments,  does 
not  indicate  appreciable  change  in  the  amplitude  of  SH 
waves  as  a function  of  the  direction  of  propagation. 

5.  The  data  in  this  paper  have  not,  as  yet,  contributed  to 

a satisfactory  theoretical  explanation  for  the  generation 
of  SH  waves  from  explosive  sources." 

Geyer  and  Martner's  second  conclusion  might  be  understood  in  terms  of 
Taylor  instability  if  one  assumed  that  at  sharp,  high-contrast  interfaces  the 
material  strength  is  low  and  variable. 

In  all  the  above  studies  one  must,  of  course,  remember  that  there  is  a 
substantial  difference  in  scale  between  small  explosions  in  soil,  and  under- 
ground weapons  tests.  For  example,  it  might  have  been  expected  that  any 
displacement  of  an  NTS  cavity  or  collapse  cone  similar  to  that  in  Figure  4 
would  have  been  noticed. 

Evidence  for  tectonic  strain  release  is  given  by  Nuttli  (1969),  who 
showed  that  there  was  a node  in  the  S wave  radiation  patterns  of  GREELEY  and 
HALFBEAK,  as  seen  in  Figures  6a  and  6b  in  which  the  direction  of  first 
motion  chang a with  azimuth. 

W.  R.  Perret  in  several  publications  has  demonstrated  the  presence  of 
substantial  close-in  SH  radiation  from  nuclear  explosions.  Figure  7,  from 
Perret  (1972)  shows  the  instrumental  layout  for  the  GASBUGGY  gas  stimulation 
experiment.  Figures  8 through  14  give  his  recorded  velocity,  displacement 
and  reduced  displacement  potential  traces.  Only  the  radial  data  (Figures  10, 
11,  and  14)  were  reproduced  In  Perret  (1972);  the  vertical  and  transverse 
data  may  be  found  in  Perret  (1970).  We  see  from  these  data  that  the  vertical 
and  transverse  motion  at  shot  level  is  typically  1/4  the  amplitude  of  the 
radial  motion. 


In  Figure  15  (from  Perret,  1968)  we  see  the  instrumental  layout  for  the 
STERLING  explosion  which  was  set  off  in  the  cavity  created  in  a salt  dome  by 
the  SALMON  explosion.  Figures  16  and  17  give  the  radial  and  vertical  shot- 
level  recordings  and  we  see  that  at  the  closest  distance  the  SV  to  P ratio  is 
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abOUt  1/3  bUt  at  sreater  dlst—  the  ratio  is  about  1/10.  Analysis  of  a 
rais-oriented  instrument  enabled  Ferret  to  eonolude  that  the  SH  to  SV  ratio 

« very  small.  These  results  are  consistent  with  the  idea  that  the  SV  radia- 
tron  observed  resulted  from  conversion  due  to  the  non-spherical  cavity 
resulting  from  the  recrystallised  pool  of  molten  salt  from  the  SALMON 
explosion.  The  absence  of  SH  radiation  is  open  to  the  natural  explanation 

that  because  of  decoupling  there  could  be  neither  strain  release.  Taylor 
instability,  nor  cracking. 

On  the  other  hand,  Ferret's  (1567)  data  for  the  original  SALMON  shot 
(Figures  18,  19,  20,  show  SH  to  P amplitude  ratios  ranging  between  0.5  and 

TheSe  data  are  di£flcult  to  understand  in  view  of  the  teleseismic 
observations  (Teledyne  Geotech,  1964)  of  no  detectible  horizontally  polarised 
long-period  radiation  from  SALMON.  He  have  verified  this  conclusion  by 
examination  of  the  original  LRSM  film  data  at  EU-AL  and  JE-LA,  the  two 
closest  stations.  At  JE-LA,  which  is  oriented  nearly  radial  and  transverse 
to  SALMON,  the  Love  to  Rayleigh  ratio  is  less  than  0.1.  The  short-period 
records,  dominated  by  2 Hs  energy,  show  a transverse  to  radial  ratio  of  about 
0.5.  One  might  argue  that  EU-AL  and  JE-LA  were' near  nodes  in  the  long  period 
Love  wave  radiation  patterns,  while  more  distant  stations  had  insufficient 
signal  to  noise  ratios  for  reliable  estimation  of  Love  waves.  Ferret's 
observations  of  S/P  - 1/3  for  explosion  together  with  the  6-7/1  ratio  expected 
rom  a simple  double-couple  earthquake  lead  one  to  expect  about  1.3  magnitude 
units  S-LR  separation,  assuming  that  S and  P convert  to  LR  equally  effici- 
ently  at  the  free  surface. 

Joksoa  and  Kehrer  (1972),  Thomson  et  al.  (1969),  Lambert  et  al.  (1972) 
Toksoz  et  al.  (1971),  ToksSz  et  al.  (1965),  Brune  and  Pomeroy  (1963)  and 
others  have  generally  interpreted  the  Love  wave  radiation  from  underground 
nuclear  tests  as  being  evidence  for  earthquake  triggering,  or  as  evidence  of 
release  of  prestress  due  to  the  bomb  cavity  or  due  to  cracks  controlled  by 
the  prestress.  It  seems  clear,  however,  in  view  of  the  results  of  Kisslinger 

et  al.  (1961)  and  Geyer  and  Mariner  (1969)  quoted  above,  that  there  are 
alternative  explanations  for  the  data. 

„F°r  "™P1C'  P-rhai'S  the  “st  convincing  evidence  for  strain  release  in 
toksoz  et  al.  (1972)  is  the  observed  reversal  of  phase  for  Rayleigh  waves 
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ho  tween  explosions  with  nniii  I 1 and  large  K values.  (Toksoz  et  al.,  and 
earlier  workers,  del  I no  !■'  as  Lhe  fraction  ol  1 aw- 1 requency  d I sp  lacement 


amplitude  due  to  sLratn  release.)  however,  as  seen  In  figure  5,  similar 
phenomena  were  I mind  by  Kissllnger  et  ai.  lor  Lhe  horizontally  polarized 


surl ace  radiation  I rom  their  explosions  In  soli  and  these  reversals  could  not 


have  been  caused  by  strain  release. 


As  a summary  ol  this  literature  review  on  the  origin  ol  SII  radiation 
I rom  underground  explosions  we  may  say  that  while  strain  release  Is  a tenable 
hypothesis  which  can  explain  anomalies  in  radiation  from  nuclear  explosions 
at  N1S;  the  existence  of  SII  radiation  from  smaller  explosions,  which  cannot 
be  explained  by  this  mechanism,  leaves  the  hypothesis  open  to  doubt. 
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DETECTION 

Table  I is  a list  of  56  events  for  which  we  attempted  to  detect  long- 
period  S and  P waves  at  LASA.  Table  II  gives  the  criteria  which  define  the 
population.  There  were  approximately  10  events  which  fit  the  criteria  but 
which  could  not  be  examined  because  of  unrecoverable  data.  Table  III  defines 
the  tour-dimensional  "window"  used  for  the  detection  study.  If  in  one  of 
the  three  128-second  overlapping  detection  windows  surrounding  rhe  expected 
arrival  time  the  program  FKCOMB  (Smart,  1972)  printed  out  a detection  meeting 
these  criteria,  then  the  event  was  declared  detected. 

The  false  alarm  rate  was  evaluated  by  looking  for  detections  in  the  other 
three  quadrants  for  those  time  windows  in  which  there  was  no  detection  in  the 
quadrant  from  which  the  arrival  was  expected.  An  S or  P in  the  FA  column  in 
Table  I indicates  the  phase  for  which  a false  alarm  was  obtained.  Some  false 
alarms  were  detected  in  runs  on  data  not  presented'here. 

Out  of  156  P time  windows  examined,  there  were  3 false  alarms.  Each 
time  window  had  three  possible  quadrants;  however,  a false  alarm  in  any  one 
of  three  time  windows  constitutes  a detection;  therefore  the  P false  alarm 
probability  per  event  is  given  by  (3/156)  x (1/3)  x 3 = 0.02.  Similarly, 
there  were  10  false  alarms  out  of  100  S wave  windows  examined  for  a false 
alarm  probability  of  0.1,  five  times  the  P-wave  value. 

A full  discussion  of  program  FKCOMB  may  be  found  in  Smart  (1972)  and 
Smart  and  Flinn  (1971);  however,  we  will  give  a brief  discussion  of  the 
method  here.  The  program  forms  the  full  three-dimensional  f-k  spectrum  and 
searches  for  three-dimensional  maxima.  It  prints  these  out  with  the 
appropriate  estimated  parameters.  The  program  has  an  option  to  strip  away 
the  energy  resulting  from  a large  interfering  signal  in  order  to  look  at  a 
smaller  one;  however,  all  our  statistics  are  without  stripping,  since  in 
running  about  20  events  with  stripping  we  did  not  pick  up  any  additional 
detections . 

To  conserve  computation  time  we  attempted  to  detect  SH  waves  only,  since 
we  thought  at  the  beginning  of  this  study  that  SH  waves  would  be  most  powerful 
for  discrimination.  An  analysis  of  long-period  SH  and  SV  waves  at  N0RSAR 
beamed  for  about  50  Asian  earthquakes  yielded  a median  SH  to  SV  amplitude 
ratio  of  0.7.  This  would  imply  that  SH  thresholds  achieved  in  this  report 
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TABLE  I 

Events  Analyzed 


Detection 


Date 

Origin  Time 

Lat(N) 

Long(E) 

Depth 

“b 

P 

S 

FA 

2 

June 

19  70 

23:33:30.2 

45.7 

150  9 

20 

5,4 

Y 

Y 

10 

June 

1970 

16:17:48.7 

44.9 

149.5 

57 

5.7 

Y 

Y 

12 

June 

1970 

03:05:20.7 

44.9 

148.9 

33 

4.9 

N 

N 

12 

June 

1970 

12:47:31.4 

43.5 

147.0 

40 

4.9 

N 

BT 

13 

June 

1970 

11:48:46.3 

44.6 

148.3 

50 

5.0 

N 

N 

22 

June 

1970 

21:33:32.6 

43.5 

147.6 

33 

5.6 

Y 

Y 

18 

Dec 

1970 

22:28:54.1 

50.2 

156.7 

61 

5.1 

Y 

Y 

21 

Dec 

1970 

10:52:54.9 

43.8 

151.0 

38 

5.1 

N 

Y 

27 

Dec 

1970 

20:44:48.9 

44.9 

150.7 

48 

5.2 

N 

Y 

6 

Feb 

1971 

10:43:05.5 

47.0 

154.1 

60 

5.6 

Y 

Y 

17 

Feb 

1971 

00:43:45.5 

46.8 

152.8 

47 

4.9 

N 

Y 

27 

Feb 

19  71 

02:26:41.8 

45.6 

150.8 

33 

5.1 

N 

N 

18 

Mar 

1971 

03:13:15.4 

49.2 

156.3 

33 

5.1 

N 

Y 

24 

Mar 

1971 

18:25:17.9 

45.0 

148.9 

56 

5.0 

N 

N 

S 

11 

June 

1971 

10:48:47.9 

49.9 

156.5 

55 

5.2 

Y 

Y 

17 

June 

1971 

09:32:0o.4 

44.4 

148.9 

33 

4.9 

N 

N 

S 

9 

July 

1971 

06:39:47.0 

43.4 

147.7 

49 

5.1 

N 

N 

s 

9 

July 

1971 

08:06:35.2 

43.4 

147.7 

54 

4.9 

N 

N 

9 

July 

1971 

16:44:15.8 

43.4 

147.7 

46 

4.9 

N 

Y 

1 

Aug 

19  71 

02:06:06.6 

50.4 

156.8 

20 

5.6 

Y 

Y 

19 

Aug 

1971 

22:15:37.7 

49,3 

155.4 

33 

6.0 

Y 

Y 

23 

Aug 

1971 

21:55:17.7 

45.6 

151.0 

34 

5.7  ■ 

Y 

Y 

p 

1 

Sept 

19  71 

10:54:04.6 

48.4 

154.9 

50 

4.9 

N 

N 

8 

Sept 

1971 

07:25:14.5 

37.2 

141.3 

56 

5.5 

Y 

Y 

9 

Sept 

1971 

10:35:23.2 

46.0 

150.0 

47 

5.1 

N 

Y 

9 

Sept 

1971 

23:01:06.8 

44.4 

150.9 

7 

6.0 

Y 

Y 

14 

Sept 

1971 

06:56:31.8 

42.3 

144,7 

39 

4.9 

N 

N 

s 

15 

Sept 

1971 

14:55:05.3 

39.1 

143.4 

17 

5.8 

Y 

Y 




TABLE  I (Continued) 
Events  Analyzed 


Origin  Time  Lat(N) 


Long(E)  Depth 


Detection 
P S FA 


TABLE  II 

Criteria  for  Selecting  Events 
from  the  U3GS  Earthquake  List 

h <_  100  km 
nib  — ^*9 

June  1970  GR  221  and  222  (Kurils) 

Dec  1971  - 30  March  1971 

June  1971  - August  1971 

nib  ^ 4. 8 

Sept  If  71  - Feb  1972  GR  221,  222,  224,  228,  229 

(Kurils  and  Japan) 
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TABLE  III 


Criteria  for  Detection  of 
Long-Period  P and  S at 
LASA  from  Kurils  and  Japan 
using  Program  FKCOMB 


Veloc ity 


Azimuth 


«mld  be  lowered  by  about  .15  magnitude  units  had  FKCOMB  been  called  upon  to 

make  the  appropriate  rotations  for  SV.  Inspection  of  these  same  data  shows 

that  SV  is  more  easily  detected  on  the  radial  than  on  the  vertical  component 

reflecting  the  fact  that  SV  motion  has  a larger  radial  than  vertical  compo-  ’ 
nent. 

The  P and  S columns  of  Table  I show  whether  or  not  the  P or  S wave  was 
detected.  Figures  21  and  22  are  computer-fitted  cumulative  normal  curves  for 
these  data  giving  the  probability  of  detecting  P or  S waves  as  a function  of 
mg.  The  90 7.  point  for  P and  S waves  are  5 6 and  5.2  respectively.  The 
50%  points  are  5.25  and  5.0,  respectively. 

We  attempted  to  evaluate  the  efficacy  of  FKCOMB  as  compared  to  beamforming 
by  beaming  P waves  from  all  magnitude  5.2  and  5.3  events.  We  also  beamed  all 
the  S wave  FKCOMB  non-detections  for  events  with  > 5.0.  These  events  are 
indicated  by  underlining  in  the  P or  S detection  columns  of  Table  I,  and  in 
every  case  if  FKCOMB  detected  or  did  not  detect  the  signal  we  reached  the 
same  conclusion  upon  viewing  the  plotted  beam.  In  general  there  seems  to  be 
a close  equivalence  between  the  power  of  the  two  techniques. 

We  have,  however,  found  other  examples  in  which  FKCOMB  failed  to  detect 
an  S wave  which  was  apparent  on  the  beam.  Figure  23  shows  the  SH  beam  for  an 
earthquake  in  Kamchatka.  The  time  windows  used  in  an  attempt  at  detection  by 

FKCOMB  are  indicated  on  the  plot.  No  suitable  detection  was  declared  on  the 
FKCOMB  printout. 

An  even  more  interesting  example  of  the  failure  of  FKCOMB  to  detect  is 
given  in  Figure  24.  FKCOMB  failed  to  detect  the  CANNIKIN  SH  wave  which  is 
quite  clear  on  the  beam.  The  reason  for  this  is,  as  we  shall  see  in  Figure  33, 
that  the  SH  component  is  quite  different  on  different  channels  at  LASA. 

FKCOMB  detections  are  determined  by  an  F statistic  which  is  large  only  if  the 
signal  is  nearly  the  same  on  every  channel. 

One  interpretation  of  this  result  is  that  the  true  SH  emitted  from  the 
CANNIKIN  source  was  negligible,  and  that  the  transversely  polarised  signals 

at  each  LASA  sensor  consisted  only  of  SV  waves  which  had  been  refracted  through 
angles  which  varied  from  site  to  site. 
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Comparison  of  radial  and  transverse  components  for  the  LR  and  LQ  waves 
fiom  CANNIKIN  show  that  the  calibrations  and  orientations  of  the  individual 
sites  were  accurate.  The  calibrations  closest  in  time  to  CANNIKIN  both 
before  and  after  the  event  were  reviewed  and  found  to  be  in  agreement  with 
each  other  to  within  5%. 


> 


DISTANCE-AMPLITUDE  RELATIONS 

As  discussed  in  the  Introduction,  the  best  previously  available  distance- 
amplitude  relation  was  presumably  that  of  Gutenberg  and  Richter  (1956),  and 
it  was  not  clear  how  accurate  it  would  be  for  use  with  modern  instrumentation 
with  different  passbands. 

We  have  therefore  developed  a preliminary  distance-amplitude  relation 
for  long-period  S from  the  data  base  developed  by  von  Seggern  (1974).  The 
data  base  consisted  of  all  long-period  S waves  measured  at  the  operational 
Long  Period  Experiment  (LPE)  stations  for  all  January-Apr il  1972  events  for 
which  LR  amplitudes  with  periods  between  17  and  23  seconds  were  measured  at 
3 or  mor w LPE  stations,  and  for  which  the  depth  was  less  than  60  km.  For 

each  such  event  an  average  value  was  calculated  using  the  formulas: 

Mg  = log  (A/T)  + 1.16  log  A + .74  A < 15° 

Mg  = log  (A/T)  + 1.66  log  A - 0.18  A > 15° 

where  A is  peak-to-peak  amplitude  in  millimicrons.  These  formulas  are  due 
to  von  Seggern  (1970)  and  to  Gutenberg  (1954a)  respectively. 

The  event  average  values  as  subtracted  from  each  log(A/T)  value  for 

S from  that  event,  and  the  results  plotted  in  2°  intervals  as  a function  of 
distance,  as  seen  in  Figure  25.  In  each  2°  interval  the  location  of  the 
mean  was  determined,  and  a smooth  curve  was  then  drawn  by  hand.  The  tabulated 

values  of  this  curve  are  given  in  Table  IV.  Use  of  this  distance-amplitude 

formula  must  be  guided  by  the  realization  that  the  seismic  analysts  picked 
the  largest  amplitude  on  either  component,  and  that  they  made  no  attempt  to 
exclude  either  late  branches  of  the  S-wave  or  PL-type  motion. 

On  Figure  25  we  have  plotted  a few  B values  from  Gutenberg  and  Richter's 

(1965)  surface  focus  distance-amplitude  relation  (see  Figure  1).  We  see  that 

the  agreement  is  probably  within  experimental  error  except  around  20°  where 
Gutenberg  and  Richter's  data  seem  to  be  too  low.  Also  in  Figure  5 we  have 
plotted  a few  points  from  the  smooth  curve  drawn  through  Nuttli's  (1969)  data 
in  Figure  6.  These  are  normalized  at  20°  and  v;e  see  that  the  agreement  is 
excellent  despite  the  fact  that  Nuttli  measured  only  the  first  peak  and  only 
for  SH. 


A 

B 

GB 

A 

B 

10 

2.61 

54 

3.50 

12 

2.50 

56 

3.58 

14 

2.48 

58 

3.71 

16 

2.45 

60 

3.75 

18 

2.48 

62 

3.76 

20 

2.50 

3.0 

64 

3.77 

22 

2.57 

66 

3.78 

24 

2.75 

68 

3.78 

26 

3.00 

70 

3.78 

28 

3.27 

72 

3.78 

30 

3.36 

3.35 

74 

3.78 

32 

3.40 

76 

3.78 

34 

3.43 

78 

3.78 

36 

3.44 

80 

3.79 

38 

3.44 

82 

3.80 

40 

3.45 

3.7 

84 

3.81 

42 

3.45 

88 

3.82 

44 

3.45 

90 

3.86 

46 

3.45 

92 

3.90 

48 

3.45 

94 

3.95 

50 

3.45 

3.6 

96 

4.00 

52 

3.46 

98 

4.10 

100 

4.15 

Jan-Apr;  LPE  Data;  h < 60  km. 

LR  detected  at  3 or  more  LPE  stations, 
GB  is  Gutenberg  B factor. 
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As  an  additional  check  on  the  accuracy  of  our  distance- amplitude  curve, 
we  have  plotted  in  Figure  26  the  points  on  Figure  25  for  those  events  for 
which  LR  was  measured  at  7 or  more  of  the  8 LPE  stations.  By  thus  restricting 
the  data  base  to  large  events,  we  hope  to  detect  any  bias  due  to  preferential 
detection  of  large-amplitude  S waves.  By  comparison  with  the  original  curve, 
also  plotted  in  Figure  26  we  see  that  it  does  not  appear  to  be  biased  with 


respect  to  the  selected  data. 


By  use  of  this  distance-amplitude  relation  we  may  attempt  to  reconcile 
our  S wave  thresholds  with  those  found  by  Strauss  (1973).  In  Figure  27  we 
see  his  results,  together  with  a fitted  cumulative  normal  curve  suggesting 


that  the  50  and  90%  points  are  m^  4.2  and  4.7  respectively.  The  range  of 
distances  from  ALPA  to  events  in  Kamchatka  and  the  Kurils  is  27°-42°,  over 


which  our  B factor  varies  from  3.13  to  3.45,  a range  of  0.32  magnitude  units. 


Let  us  assume  an  average  B factor  of  3.40,  keeping  in  mind  the  sharp  increase 
in  the  B factor  right  at  27°  and  the  greater  seismicity  in  the  Kuril  Islands. 

The  range  of  distances  from  LASA  to  events  in  the  Kurils  and  Japan  is  56°  to 
84°  with  a range  of  B factors  between  3.58  and  3.81.  Over  most  of  this  range 
the  B factor  is  3.78.  Let  us  select  this  as  an  average  B factor.  The  dif- 
ference in  average  B factors  is  0.38  magnitude  units.  Since  Mg  is  approximately 
proportional  to  the  first  power  of  m^  in  Seismic  Region  19  (von  Seggern,  1974) 
this  difference  in  long  period  magnitude  should  translate  into  an  approximately 
equal  difference  in  short-period  thresholds „ 


From  Figures  22  and  27  we  see  that  the  observed  difference  in  50% 


thresholds  is  0,78  m^  units,  leaving  a difference  of  0.4  magnitude  units 
unexplained. 


This  difference  may  be  accounted  for  by  the  differing  noise  levels  at 
the  two  arrays.  From  Strauss  (1973,  Figures  1,  2,  3)  we  can  estimate  the  rms 
18-40  second  median  single  sensor  LPZ  noise  level  to  be  10  mp.  Dividing  by 
the  square  root  of  19  to  obtain  the  beam  rms  amplitudes  we  obtain  2,3  mp 
(from  Clark  et  al.  1972,  Table  V,  one  may  obtain  the  median  ALPA  LPZ  beam 
rms  value  as  7.5  mp,  but  their  bandpass  limits  were  15-50  seconds,  and  the 
data  are  from  1970,  when  there  was  a severe  long-period  system  noise  problem 
at  ALPA).  The  noise  level  at  LASA  may  be  obtained  from  Farnham  (1968), 


V' 
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Table  II,  and  Masse  et  al.  (1970,  Table  III).  The  median  rms  LPZ  15-50  second 
value  from  these  studies  is  4.5  my.  Thus  we  see  that  the  beam  noise  level  at 
ALPS  Is  approximately  4.572.3  - 2.0  times  lower  than  at  I .ASA.  Assuming  that 
the  horizontal  noise  levels  are  proportional  to  the  vertical  noise  levels, 
this  could  account  for  0.3  of  the  0.4  magnitude  units  discrepancy.  The 
discrepancy  could  be  further  reduced  to  0.05  magnitude  units  by  allowing  for 
the  SH/SV  ratio  of  -0.15  magnitude  units  as  discussed  in  the  previous  section. 
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MEASUREMENTS  OF  S FROM  EXPLOSIONS,  AND  DISCRIMINATION 
WITH  RESPECT  TO  EARTHQUAKES 


fl 
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In  measuring  S waves  from  NTS  explosions  we  decided  to  concentrate  on 
those  LRSM  stations  oriented  radially  and  transverse  to  the  test  site  in 
order  to  eliminate  the  possibility  that  inaccurate  calibrations  could  intro- 
duce an  apparent  SH  component  which  was  in  fact  rotated  SV  motion.  The 
further  requirement  that  the  stations  be  in  operation  for  most  of  the  events 
of  interest  reduced  the  list  of  stations  to  RKON,  WH2YK,  !INME,  and  NPNT. 

For  MILROW  and  CANNIKIN  in  this  report,  we  have  used  some  rotated  data,  ana 
for  the  CaNNIKIN  event  several  other  LRSM  stations  were  specially  deployed 
radial  and  transverse  to  the  event . 

For  the  sake  of  concreteness,  we  have  included  in  this  report  several 
illustrations  of  S waves  from  explosions  Figure  28  gives  the  rotated 
vertical,  radial,  and  transverse  waveforms  from  CANNIKIN  as  observed 
at  the  LASA  AO  subarray,  A distinct  SH  component  about  1/4  the  amplitude  of 
SV  can  be  seen.  Figures  29,  30,  and  31  give  a more  detailed  look  at  the 
vertical,  radial,  and  transverse  components  of  the  CANNIKIN  P wave.  Note 
the  very  small  upward  first  motion,  which  might  well  be  attributed  to  the 
effect  of  pP.  (It  is  possible  that  the  P waveforms  are  distorted  by  non- 
linearities  in  the  amplifier,  which  will  clip  at  30  mv  peak-to-peak  input 
(Gudzin  and  Hennen,  1967  and  Gudzin,  personal  communication).  The  seismometer 
output  is  4.5  mv/u  at  1 Hz,  and  the  CANNIKIN  short-period  P at  LASA  is  about 
6u  peak-to-peak  at  .9  Hz,  which  yields  24  mv  peak-to-peak.)  Note  that  the 
transverse  component  is  equal  in  amplitude  to  the  radial  although  both  are 
only  about  0.1  times  the  vertical.  Figures  32  and  33  give  the  radial  and 
transverse  components  of  the  CANNIKIN  S wave  at  all  the  LASA  subarrays.  Note 
that  the  amplitude  for  SV  ranges  over  a factor  of  2460/1504  = 1,6.  A review 
of  the  instrument  calibrations  shows  that  the  variation  must  be  real.  The 
variation  is  even  greater  for  SH,  on  the  order  of  1064/341  = 3,1,  and  it  is 
clear  that  the  actual  waveform  varies  in  shape  from  instrument  to  instrument. 
Figures  34  and  35  show  the  SV  and  SH  waves  at  LASA  for  MILROW,  We  see  that 
the  beam  SH/SV  ratio  is  much  greater  for  MILROW  than  for  CANNIKIN.  Figures  36- 39a 
show  tracings  of  the  best  signal  - to-noise  ratio  data  available  at  RKON,  WH2YK, 
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HNME,  and  NPNT . In  Figures  38  and  39a  for  HNME  and  NPNT  we  see  that  following 
an  initial  motion  on  R and  T,  the  motion  grows  larger  and  appears  on  Z»  This 
motion  may  be  interpreted  either  as  PL  motion  (Oliver,  1961;  Helmberger  and 
Engen,  1974);  or  as  late  arrivals  due  to  triplication  of  the  travel-time 
curve  (Nuttli , 1959,  1964;  Nuttli  and  Whitmore,  1962;  Haskell,  1962).  Note 
that  at  HNME  this  late-arriving  energy  also  shows  up  clearly  on  the  transverse 
component,  suggesting  that  if  it  is  due  to  PL-type  conversion,  then  it  must 
be  converted  from  the  incoming  SH  motion,  a possibility  suggested  by  Oliver 
(1961). 

For  routine  discrimination  purposes  it  seems  clear,  considering  the 
controversy  over  the  nature  of  this  energy,  that  we  must  rely  on  the  analyst 
only  to  pick  the  maximum  amplitude  in  the  first  few  cycles  of  an  apparent 
phase.  This  is  the  approach  taken  in  the  present  report. 

In  Figure  40  we  show  tracings  of  S waves  at  RK.0N  for  PILEDRIVER. 
PILEDRIVER  is  of  special  interest  because  it  has  the  highest  F value  reported 
by  Toksoz  and  Kehrer  (1972). 

In  selecting  events  to  analyze  we  decided  to  restrict  ourselves  to 
explosions  of  about  a megaton  or  larger  in  order  to  avoid  signal-to-noise 
ratio  problems  and  to  ensure  that  we  fully  sampled  the  range  of  S wave  ampli- 
tudes emitted  from  the  explosions  The  resulting  set  of  events  is  described 
in  Table  V.  FAULTLESS  is  included,  even  though  the  yield  is  listed  by 
Springer  and  Kinnaman  (1971)  as  intermediate,  because  it  was  exploded  at  a 
depth  greater  than  3000  feet,  and  because  Its  Mg  value  is  larger  than  all 

other  events  in  the  list  except  BOXCAR.  The  M values  for  FAULTLESS.  HANDLY , 

s 

and  JORUM  were  determined  by  us,  using  LR  readings  at  NPNT,  RKON,  PGBo,  and 
HNME  together  with  the  distance  formula  and  station  corrections  given  by 
von  Seggern  (1973),  and  with  a further  correction  as  described  below. 

Von  Seggern's  formulas  were  derived  using  continental  United  States  and 

Canada  LRSM  data  for  NTS  explosions  His  station  corrections  and  distance- 

amplitude  relations  ensured  that  approximately  the  same  M value  would  result 

s 

from  measurement  of  the  maximum  film  amplitude  at  any  station.  The  overall 


, 
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level  was  set  to  yield  the  same  Mg  as  would  Gutenberg’s  formula  at  distances 
greater  than  20°.  This  was  accomplished  by  adjusting  the  overall  magnitude 
level  tu  be  in  agreement  with  that  of  Lambert  ind  Alexander  (1971).  However, 
while  almost  all  of  the  LRSM  data  even  at  distances  of  40°  in  North  America  ’ 
have  maxima  at  the  Airy  phase,  such  pure  continental  data  was  rare  in 
Gutenberg  s original  data  set  (Gutenberg,  1945a)  and  in  von  Seggern's  (1974) 
data  set,  from  which  Mg  values  for  normalization  of  our  S wave  distance- 
amplitude  relation  were  taken.  An  example  of  such  an  Airy  phase  may  be  seen 
on  the  bottom  trace  in  Figure  39b.  Note  the  strikingly  different  waveforms 

for  events  of  the  same  yield  recorded  at  the  same  station  at  approximately 
equal  distances. 

The  existence  of  this  Airy  phase  will  cause  earthquakes  or  explosions 
near  NTS,  recorded  in  North  America,  to  have  too  large  an  M . This  effect 
was  recognized  by  Marshall  and  Basham  (1972).  However,  their  correction 
factors,  derived  from  WWSSN  data,  do  not  seem  to  be  large  enough  to  account 
for  the  discrepancies  illustrated  in  Figure  39b.  This  may  be  due  to  the  fact 
that  because  of  the  narrower  response  of  the  LRSM  system,  the  Airy  phase 
maximum  occurs  at  a longer  period  than  on  the  WWSSN  system.  Another  explana- 
tion would  be  that  the  dispersion  curves  chosen  by  Marshall  and  Basham  for 
North  America  are  not  truly  typical. 

We  choose  instead  to  make  our  own  correction  by  measuring  at  a period 

away  from  the  Airy  phase,  as  did  Gutenberg,  and  using  the  resulting  difference 

in  magnitude  as  a path  correction  for  magnitudes  computed  from  measurements 

made  on  the  Airy  phase.  In  reference  to  Figure  39b  we  plotted  the  group 

velocity  as  a function  of  period,  and  found  that  a period  of  23  seconds 

is  well  away  from  the  group  velocity  minimum  (dU/'dT  = 0).  We  have  therefore 

measured  the  amplitude  in  Figures  39b  at  23  seconds  period.  The  result  is 

0.53  units  smaller  than  the  magnitude  resulting  from  measurements  on  the 

Airy  phase  "<a:  imum.  The  other  megaton-sized  NTS  explosions  recorded  at  NPNT 

have  very  similar  waveforms,  and  measurements  of  the  corresponding  cycles 

yields  a mean  AM  of  0.53. 
s 

Work  by  von  Seggern  (1974)  and  other  authors  has  established  that  the 
proper  worldwide  distance  amplitude  relationship  for  LR  in  the  absence  of 
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dominant  Airy  phases  is  logA  ~ (0,9-1.0)logA  for  5°  < A < 150°.  Von  Seggern's 
1973  formula,  with  logA  ~ 0.87  logA  is  used  in  this  report  and  is  consistant 
with  the  above  values.  Thus  we  may  interpret  von  Seggern’s  station  corrections 
(<V  as  being  predominantly  due  to  effects  peculiar  to  the  Airy  phase,  and 
we  would  expect  that  they  would  not  be  satisfactory  corrections  for  measurements 
at  other  locations  on  the  waveform.  Thus,  the  C.  corrections  (0.01  for  NPNT) 
should  not  be  applied  directly  to  the  NPNT  measurements  at  23  seconds  discussed 
above , But  we  can  use  them  with  the  following  formulation: 

Gutenberg  Airy  ^Gutenberg  ^Airy ^ 

= MAiry  + {1°S(A/T>t=73  “ [1°S (A/T)Airy  “ Ci]} 

= MAiry  “ {[l0g(A/T>Airy  ~ log (A/T)^]  " C.}. 

Thus  we  shall  subtract  A = 0.52  = 0.53  - 0.01  magnitude  units  from  all 
Mg  values  determined  using  the  method  of  von  Seggern  (1973).  As  a check  on 
this  result,  we  analyzed  JORUM  at  HNME  where  C..  = -0.20  and  found 
A = 0.37  - (-0.20)  = 0.57.  The  result  at  NPNT  is  preferable,  however,  since 
the  waveform  there  shows  the  clearest  dispersion.  Ideally,  we  should  deter- 
mine A for  all  stations  used  and  average  them;  but  this  is  not  possible  for 
the  other  LRSM  stations  because  their  range  from  NTS  is  too  short  for  an 
adequately  dispersed  wave  to  clearly  show  20  second  cycles. 

This  result  illustrates  a possible  general  approach  to  the  measurement 
of  Mg  and  was  analyzed  in  some  detail  by  von  Seggern  (1971).  We  suggest 
that  a graph  of  group-velocity  vs  period  be  constructed  for  each  epicentral 
region-station  pair.  A period  near  20  seconds  well  away  from  any  Airy  phases 
should  be  selected,  and  the  amplitude  corrected  using  the  formula  from 

Marshall  and  Basham  (1972)  appropriate  for  stationary  phases: 

-3/2  -1/2 

A ~ UT  (dU/dT)  . The  difference  between  the  corrected  amplitude  and 
the  maximum  amplitude  on  the  record  can  be  used  as  a path  correction  for 
measurements  made  at  the  maximum  amplitude  on  future  records.  In  addition, 
the  corrected  amplitudes  could  be  plotted  as  a function  of  distance  to 
determine  an  improved  distance-amplitude  relation. 


In  Table  VI  wt  give  the  amplitude  and  period  measurements  for  S together 
with  the  Mgy  and  M magnitudes  calculated  using  the  distance-amplitude 
formula  in  Table  IV  and  Figure  25. 


Finally,  in  Figure  4la  we  present  a plot  of  M - M versus  M - M 

oil  S j V S 


for  the  data  measured  in  this  report.  We  expect  the  earthquakes  fo  cluster 
0.1  - 0.2  magnitude  units  diagonally  below  (0,0)  since  we  defined  the 


distance-correction  terms  for  S-wave  Mg  using  measurements  taken  from  the 


component  with  the  larger  S amplitude. 


We  note  that  the  center  of  gravity  for  the  NTS  explosion  points  lies 
slightly  below  the  diagonal  line  near  the  point  (-.6,  -.9),  showing  good 
separation  from  the  mean  earthquake  population,  While  GREELEY,  with  an  F 
value  of  1.6,  is  above  the  line  for  NPNT  and  HNME,  the  absolute  value  of  the 
S wave  is  still  small  compared  to  the  Rayleigh  wave.  The  same  remark  is  true 
for  the  single  reading  at  RKON  of  an  S wave  from  PILEDRIVER,  which  has  an  F 


value  of  3.20. 


There  are  several  possible  explanations  for  the  separation  between 
earthquakes  and  explosions  other  than  energy  variation  at  the  source. 


o A low  Q region  beneath  nuclear  test  sites  absorbs  the  long-period 
S waves  but  is  deep  enough  not  to  attenuate  the  surface  waves. 


o Due  to  a low-velocity  surface  layer,  the  incidence  angle  (with 
respect  to  the  vertical)  for  S-waves  observed  at  A = 20°-40°  is  substantially 
less  than  the  value  of  30°  given  by  Chandra  (1972).  Then  if  the  strain 
release  fault  mechanism  is  vertical  strike-slip,  as  suggested  by  Toksoz 
and  Kehrer  (1972),  the  shear  radiation  would  be  close  to  the  downward  null 


axis  of  the  source  mechanism. 


o The  strain  release  fault  may  in  fact  not  be  vertical  strike-slip, 
but  may  dip  to  the  northeast,  creating  a node  in  the  tcleseismic  radiation 
pattern  for  the  signals  at  almost  all  of  our  observing  stations. 


We  will  now  discuss  each  of  these  possibilities  in  more  detail.  If  a 
low-Q  layer  existed,  one  would  expect  that  teleseismic  long-period  S arriving 
in  the  Basin  and  Range  would  have  lower  amplitudes  than  waves  from  the  same 
event  arriving  on  the  East  Coast.  However,  examination  of  the  long-period 
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s waves  for  the  same  events  and  stations  used  by  Der,  Masse,  and  Gurski 
a974)  to  analyze  short-period  S,  shows  less  than  a 10%  variation  between 

stations.  This  would  seem  to  rule  out  a low-Q  zone  as  a source  of  low  long- 
period  S waves. 

While  a low-velocity  surface  layer  might  result  in  a steeper  emergence 
angle  for  short-period  S waves,  one  would  not  expect  such  an  effect  for  20- 
second  waves  with  wavelengths  on  the  order  of  60-100  km.  Therefore  the  angles 
given  by  Chandra  are  appropriate. 

A dip  of  the  strain-release  fault-plane  to  the  northeast  cannot  be  ruled 
out  with  the  present  data.  It  is,  however,  a rather  adhoc  explanation.  Also, 

sS  would  not  be  near  a node  in  this  case  and  should  contribute  heavily  to 
the  signal. 

Not  only  is  the  overall  level  of  the  S wave  excitation  for  explosions 
lower  than  that  to  be  expected  for  earthquakes,  but  also  the  relative  ampli- 
tudes of  S waves  recorded  at  RKON  from  NTS  shots  cannot  be  predicted  by  use 
of  Toksoz  and  Kehrer's  (1972)  fault  azimuths  and  F factors. 

If  we  assume  that  M , or  a "tectonic"  MT  given  by 

s 

T p 

Ms  = 1o810(i£f  10  S> 

adequately  measures  the  strength  of  the  released  strain,  and  that  the  shear 
radiation  pattern  is  as  described  by  Ben-Menahem  et  al.  (1965)  and  Helmberger 
(1974),  then  we  may  form  the  quantity  - (M*  + log | cos  20 j ) where  0 is  the 
angle  between  the  strike  of  the  vertical  strike-slip  earthquake  as  given  by 
Toksoz  and  Kehrer  (1972)  and  the  epicenter-station  azimuth.  This  quantity 
Should  be  nearly  the  same  for  all  explosions  from  NTS  observed  at  a single 
station.  Effects  due  to  distance  may  be  neglected,  since  the  distance  from 
NTS  to  any  single  station  is  nearly  constant.  In  Table  VII  we  see  that  at 
RKON  this  quantity  is  no  more  constant  using  M^  + log | cos  20 | as  the 
normalizing  factor  than  using  Ms  itself,  and  this  suggests  that  the  SH  body 
wave  radiation  does  not  fit  naturally  into  the  double-couple  vertical 
strike-slip  source  model  suggested  by  Toksoz  and  Kehrer.  Considering  the 

expected  variance  of  the  data,  however,  we  do  not  have  a very  strong  test  of 
the  hypothesis. 
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We  mention  that  one  possibly  large  error  in  Table  VII  is  'that  our 
original  measurement  of  for  PILEDRIVER  may  have  been  seriously  too  small 
because  of  the  cancellation  of  the  direct  explosion  LR  by  the  strain  release 

ii 

(Toksoz  and  Kehrer,  1972).  However,  if  M were  made  larger  for  PILEDRIVER 

a RK  T 

the  variance  about  the  mean  of  the  numbers  in  the  M - M - log  I cos  0|  column 

a RK 

of  Table  VII  would  increase  while  decreasing  in  the  column  M - M . 

SH  s 

Before  a definite  conclusion  can  be  reached  that  these  data  do  not  support 

a double-couple  source  a more  careful  analysis  is  needed.  For  example,  it 

would  be  desirable  to  correct  the  LR  measurements  at  each  observing  station 

for  each  event  for  the  effects  due  to  strain  release  in  order  to  determine 
T 

a more  accurate  M . 

s 

To  gain  further  insight  into  the  problem  we  decided  to  analyze  earth- 
quakes near  NTS.  We  chose  all  events  between  1962  and  1970  on  the  USGS  seis- 
micity list  within  a circle  of  radius  2°  around  BENHAM,  with  m^  >_  5.0  which 
were  not  aftershocks  of  either  an  explosion  or  another  earthquake.  We  found 
three  events.  In  Table  VIII  are  the  source  parameters  for  two  of  them.  The 
third,  a USGS  m^  5.0  event  on  October  23,  1970,  had  long-period  LR  and  S radi- 
ation so  weak  that  it  could  not  be  analyzed.  In  Table  IX  we  see  the  calculated 
Ms’  MSH’  and  MSV  va^ues«  anc*  t*le  results,  plotted  on  Figure  41a,  place  the 
earthquakes  on  the  border  between  the  earthquakes  and  explosions.  M values 
were  calculated  with  the  same  procedures  used  for  NTS  explosions  including 
the  -0.52  dispersion  correction. 

As  a parallel  effort  we  analyzed  two  of  the  four  NOS  events  with  NOS  M 

s 

values  between  5.1  and  5.3  between  1962  and  1970  which  lay  within  2°  of 
MILROW,  and  which  are  closest  in  time  to  MILROW.  These  events  are  listed  in 
Table  VIII.  The  magnitude  restriction  was  imposed  in  order  to  ensure  good 
signal-to-noise  ratios  without  clipping. 

The  measurements  and  magnitude  estimates  are  given  in  Table  IX,  and  the 

results  are  plotted  in  Figure  41a  for  those  stations  for  which  components 

were  oriented  radial  and  transverse  to  Amchitka.  We  see  that  these  earthquakes 

are  typical  of  the  world-wide  population  (Gutenberg's  distance-amplitude 

formula  was  used  to  calculate  M from  LR) . 
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Except  for  the  signal  recorded  at  SJTX,  we  see  that  CANNIKIN  amplitudes 
are  in  relatively  good  agreement  with  NTS  events.  The  signal  recorded  at 
SJTX  has  a signal-to-noise  ratio  of  almost  10,  and  there  can  be  little  doubt 
as  to  its  reality.  MILROW  has  a number  of  relatively  large  S waves,  and  at 
most  of  the  stations  the  signal-to-noise  ratio  was  large  enough  to  make  mea- 
surement unambiguous.  However  the  most  anomalous  measurement,  at  KNUT,  was 
at  a low  signal-to-noise  ratio;  and  so  may  well  be  biased  toward  too  high  an 
amplitude . 

The  rather  striking  difference  between  MILROW  and  CANNIKIN  is  hard  to 
understand.  CANNIKIN  was  exploded  in  a fifty-foot  diameter  mined  cavity, 
whereas  MILROW  was  simply  lowered  down  a borehole  (N.  Pruvost,  personal 
communication).  This  suggests  that  CANNIKIN  with  its  low-amplitude  Love  waves 
(see,  e.g.,  Willis  et  al.  1972),  and  S waves  as  seen  in  this  report  may  be 
the  anomalous  event,  and  that  MILROW  may  be  the  more  normal  event  for  the 
Amchitka  test  site.  Engdahl  (1972)  has  suggested  that  there  is  very  little 
tectonic  stress  near  the  surface  of  the  crust  at  Amchitka. 


As  a preliminary  investigation  into  the  performance  of  the  S/LR  discrimi 
nation  for  Asian  explosions  we  beamed  NORSAR  on  the  large  August  28,  1972 
event  at  Novaya  Zemlya.  The  resulting  point  is  plotted  on  Figure  41a.  The 
normalizing  Mg  was  determined  by  von  Seggern  (1974)  using  LPE  data.  We  see 
that  it  fits  nicely  into  the  NTS  population.  Further  investigation  of  Asian 
explosions  and  earthquakes  is  underway  at  this  time. 


To  attempt  to  reduce  the  scatter  in  Figure  41a,  we  have  averaged  the 

individual  station  Mcu  and  M values  for  each  event.  (If  in  a series  of 

bn  b V 

M or  M measurements  for  a given  event  there  were  noise  measurements  which 
SH  SV 

yielded  low  upper-limit  magnitude  estimates,  then  the  average  magnitude  was 
replaced  by  the  median  magnitude  with  the  low  noise  measurements  included  in 
the  histogram.)  The  results  are  plotted  in  Figure  41b.  We  see  that  the  best 
line  we  can  draw  still  misclassif ies  the  MILROW  explosion  and  that  BENHAM  is 
perilously  close  to  misclassif ication . 


All  the  sample  earthquakes  fall  below  the  point  (0,0)  on  Figure  41L,  as 
would  be  expected  since  this  normalizing  point  corresponds  to  MSHEAR  determined 


fr°”  the  mSKimU°  0bserved  Shcar  on  either  component,  as  discussed  pre- 

viously. To  examine  the  discrimination  which  would  result  If  sue!  procedure 
had  been  followed  for  the  earthquakes  and  explosions  analysed  in  this  report 
we  have  computed  the  average  magnitude  for  each  event  obtained  by  averaging  ’ 
the  maximum  magnitude  component  at  each  station.  Since  SH/SV  does  not  appear 
from  figures  41a  and  41b  to  be  a good  discriminant,  little  should  be  lost  by 
this  procedure.  The  result  Is  plotted  in  Figure  41c.  He  see  that  not  only 
do  the  earthquakes  now  cluster  about  the  origin  as  they  should,  but  the 
discrimination  has  Improved  such  that  perfect  discrimination  is  possible  with 
a suitably  chosen  line.  It  must  be  admitted,  however,  that  the  populations 

are  so  close  that  a false  classification  rate  of  several  percent  must  be  ' 
expected  in  practice. 
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APPLICATION  OF  NEGATIVE  DISCRIMINANTS 


If,  as  Figure  41  suggests,  there  is  ^0.5  magnitude  difference  between 
the  amplitudes  ratio  S/LR  for  explosions  and  earthquakes,  then  application  of 
this  technique  for  positive  discrimination  (measurement  of  the  LR  and  S 
amplitudes  are  used  to  determine  the  nature  of  the  event)  for  explosions  will 
be  severely  limited.  Therefore  we  investigate  here  some  of  the  general 
requirements  for  the  application  of  negative  discriminants  in  which  the 
absence  of  detectible  S waves  is  an  indication  of  an  explosion.  The  theory 
as  developed  will  have  application  not  only  to  negative  discrimination  by 
S/LR,  but  also  by  M^m^. 

One  may  say  that  negative  discrimination  is  possible  when  no  stations 
detect  S waves  from  an  event.  We  must  note,  however,  that  in  this  case  if  S 
waves  are  not  detected  from  an  earthquake  then  there  is  a false  alarm.  Routine 
network  evaluation  programs,  e.g.  NETWORTH  (Wirth,  1970),  can  calculate  the 
probability  of  a zero-station  detection  and  can  therefore  determine  the 
threshold  magnitude  for  any  desired  false  alarm  rate.  Thus,  to  determine 
the  capability  of  negative  discrimination  it  is  only  necessary  to  compare 
the  threshold  magnitude  for  some  reasonable  false  alarm  rate  with  the  detec- 
tion threshold  required  for  positive  discrimination. 


Let  us  imagine  that  an  event  has  occurred  at  50°  N,  155°  E (Kamchatka). 
Let  us  further  suppose  that  our  25-station  long-period  network  is  similar  to 
that  of  Romney  (1971),  and  has  zero-to-peak  median  long-period  noise  levels 
of  19  mp  except  that  it  has  15  rap  noise  levels  at  LASA,  ALPA,  and  NORSAR. 

Let  us  also  assume  a standard  deviation  of  0.3  magnitude  units  in  the  signal 
and  noise  amplitudes  - and  a requirement  of  a signal- to-noise  ratio  of  1.5 
for  detection.  Then  program  NETWORTH  enables  us  to  calculate  the  probability 
of  detection  of  S waves  by  K or  more  stations  as  a function  of  magnitude. 

The  results  are  plotted  in  Figure  42. 


We  see  that  the  magnitude  for  a false  alarm  probability  of  0.01  is  almost 
exactly  equal  to  the  magnitude  for  90%  probability  of  detection  by  two  or 
more  stations.  Thus  the  M negative  discrimination  threshold  for  explosions 
is  about  equal  to  the  positive  discrimination  threshold  for  earthquakes.  The 
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positive  discrimination  threshold  for  explosions  is  0.5  Mg  units  higher  than 


for  earthquakes  because  of  the  mean  separation  of  the  two  populations  as 
seen  in  Figure  41a, b,c. 


Some  discussion  is  required  of  the  particular  probability  levels  we  have 
chosen  for  false  alarm  rate  and  probability  of  detection.  Since  there  are 


several  thousand  worldwide  Ms  - 3.4  or  larger  earthquakes  each  year;  there 


will  be  several  tens  of  false  alarms.  However,  this  technique  is  only  to  be 
regarded  as  a screening  mechanism,  and  each  of  the  events  could  be  investi- 
gated further.  In  addition,  the  chief  application  of  the  technique  will  be 
against  shot  arrays;  and  the  possible  locations  for  these  have  only  several 
hundreds  of  candidate  earthquakes  per  year. 


The  requirement  for  detection  by  2 or  more  stations  has  several  justi- 
fications : 


o The  probability  is  reduced  of  picking  an  arrival  from  a 
different  event  not  detected  by  the  short-period  network, 
o The  amplitude  of  an  isolated  detection  may  be  anomalously  large 
because  of  multipath  focussing  or  other  unknown  cause;  and  may 
result  in  an  explosion  being  classified  as  an  earthquake. 

We  should  note  that  the  differences  between  the  thresholds  for  negative  and 
positive  discrimination  would  be  substantially  reduced  if  the  requirement 
were  for  .001  probability  of  a false  alarm  and  90%  probability  of  detection 
by  1 or  more  stations.  The  conclusions  would,  of  course,  be  completely 
reversed  were  we  to  look  at  the  question  through  the  eye  of  a cautious  evader 
and  not  allow  even  .01  probability  of  detection  by  1 or  more  stations. 

In  Figure  42  we  see  that  the  90%  threshold  for  detection  by  4 or  more 
stations  is  3.62,  about  0.2  magnitude  units  higher  than  the  thresholds 
discussed  above.  The  same  relationships  between  thresholds  are  generally 
valid  at  all  points  on  the  globe,  as  can  be  seen  by  comparing  Figures  43-45, 


Figure  43  gives  contours  of  the  Mg  magnitude  for  which  there  is  a .01 


probability  of  a false  alarm.  Figures  44  and  45  give  contours  of  the  M 
magnitude  for  which  there  is  a 90%  probability  of  detection  by  2 or  4 stations 


i 


respectively.  The  star  on  each  of  these  Figures  indicates  the  epicenter  for 
which  figure  42  was  calculated.  The  triangles  indicate  the  station  locations. 
One  station  is  off  the  map  at  the  South  Pole, 


While  the  standard  deviations  of  t),3  magnitude  units  used  for  the 
foregoing  calculations  are  reasonable  for  ensembles  of  noise  amplitudes  mea- 
sured over  the  course  of  a year  (von  Seggern,  1974)  and  for  S wave  signals 
with  respect  to  their  mean  (see  Figure  25)  there  are  certain  situations  where 
0.1  magnitude  units  is  more  reasonable;  we  next  investigate  the  relationship 
between  thresholds  in  this  case. 

Such  a situation  arises  when  c-ne  considers  discrimination  of  a particular 
event  in  a well-studied  location.  In  this  case  the  rms  noise  amplitude  is 
known  exactly,  and  station-epicenter  amplitude  anomalies  may  have  been  deter- 
mined. From  von  Seggern  and  Blandford  (1974),  von  Seggern  (19  73),  and 
Blandford  and  Wirth  (1973),  it  is  then  reasonable  to  assume  a standard  devia- 
tion of  0.1  magnitude  units  for  both  noise  and  signal. 

To  look  at  such  a situation  from  another  point  of  view,  consider  the 
case  of  evaders  attempting  to  avoid  disc  Mmination  by  the  S/LR  criterion. 

We  may  presume  that  they  know'  the  rms  noise  levels  at  each  detecting  station, 
as  well  as  the  station  amplitude  anomalies.  They  could  therefore  attempt 
very  precisely  to  just  slip  under  the  detection  threshold. 

In  Figure  46  we  see  the  curves  of  0.1  standard  deviations  corresponding 
to  Figure  42.  We  note  that  the  curves  are  much  steeper  and  that  the  magnitude 
for  0.01  probability  of  a false  alarm  is  about  0.2  unit  below  the  magnitude 
for  90%  probability  of  detection  by  2 or  more  stations.  Thus  the  negative 
discrimination  threshold  for  explosions  is  0.2  magnitude  units  lower  than 
the  positive  discrimination  threshold  for  earthquakes.  The  positive  discrimi- 
nation threshold  for  explosions  is  0,7  = 0.2  + 0.5  units  higher  than  for 
earthquakes  because  of  the  mean  separation  of  the  two  populations  as  seen  in 
Figure  41.  Inspection  of  Figures  47-49  will  show  that  conclusions  drawn  from 
Figure  46  are  qualitatively  valid  at  most  geographical  locations.  Furthermore, 
while  these  figures  are  in  general  similar  to  Figures  43-45,  they  are  more 
complicated  and  show  sharper  gradients  in  threshold  magnitude.  This  is  because 
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the  sharp  maximum  in  the  distance- amplitude  relationship  of  Figure  25  has  not 
been  "smoothed  out"  by  the  large  values  of  a.  For  example,  in  Figure  47  com- 
pare the  deep  minimum  around  15°  S,  150°  E with  the  corresponding  point  in 
Figure  43,  and  compare  the  minima  in  Figure  49  with  the  corresponding  positions 
in  Figure  45. 


The  existence  of  an  isolated  station  can  lead  to  very  large  threshold 
differences  if  o = 0.1.  Note  the  difference  of  0.8  magnitude  units  between 


Figures  47  and  48  at  75  S,  180  E.  This  is  due  to  the  isolated  station  at  the 
South  Pole.  In  this  case  the  negative  threshold  for  explosions  is  1.3  magni- 
tude units  below  the  explosion  positive  threshold, 


i igures  50  57  tell  a similar  story  for  detection  of  Rayleigh  waves  for 
use  in  positive  and  negative  discrimination  by  use  of  LR  in  M :im  . In  these 
calculations  we  used  the  distance-amplitude  formula 


Mg  = log1()(A/T)  + log  A + 1.12 


where  A is  peak-to-peak  amplitude  in  my,  and  T was  taken  to  be  20  seconds. 
For  A < 10  degrees  the  B factor  was  taken  equal  to  its  value  at  10  degrees 
in  order  to  avoid  concentric  threshold  contours  around  each  station  for  the 
case  of  detection  by  one  or  more  stations.  The  main  difference  observed 
between  LR  and  S waves  is  that  for  LR  we  do  not  find  the  dramatic  spread  in 


thresholds  for  0 = 0.1  that  we  found  for  S, 


While  we  have  found  great  value  in  the  concept  of  negative  discrimination, 
it  must  be  remembered  that  the  theory  relies  on  an  extrapolation  of  the 
distribution  of  recorded  amplitudes  to  probability  levels  such  as  0.01,  0.001, 
where  the  supporting  data  bases  are  just  beginning  to  reach.  Furthermore, 
the  application  of  negative  discrimination  in  practice  requires  careful  atten- 
tion to  details  of  system  failure,  a very  careful  review  of  all  records,  and 
careful  noise  amplitude  measurements  at  selected  periods  by  experienced  ana- 
lyst.'. Finally,  the  resulting  discrimination  is  often  felt  to  be  unsatisfac- 
tory for  various  reasons  probably  closely  related  to  these  difficulties. 


Further  study  of  the  actual  operation  of  a network  using  negative  and 
positive  discriminants  is  currently  planned  by  using  the  network  simulation 
capabilities  recently  developed  by  von  Seggern  and  Blandford  (1974). 


N, 


hr- 


-40- 


SUMMARY  AND  SUGGESTIONS  FOR  FURTHER  RESEARCH 
We  have  developed  a distance- amplitude  formula  for  S waves  and  found  at 

» g°°d  °£ <i969)  and  the  d;n;i;e  out 

amplitude  relation  o£  Cutenbetg  and  Kichter  (1956).  except  near  20  out 

„ed  This  amplitude  maximum  suggests  that  stations  in 

maximum  is  more  pronounced.  Ihisampll  , . „/LE 

the  distance  range  15"-25"  would  be  very  valuable  for  application 
ratio  discriminant. 

We  find  the  50%  and  90%  thresholds  tor  detection  of  S waves  at  bASA  to 
b “ and  5.2  respectively,  using  either  beamforming  or  PKCOHB.  These 

v ll  are  cListent  with  the  values  ,2  and  9,  from  Strauss,  ,1923)  AIM 
data  when  consideration  is  taken  of  the  noise  levels  and  distance-amplitude 

relation  for  S waves. 

Application  of  our  distance-amplitude  relationship  made  it  possible  to 
1 ate  M and  M versus  M as  a discriminant,  and  the  results  s 
C1I  .J!  b,  c, "suggest  thSat  it  is  a good  discriminant  Such  a discriminant 
should  worir  as  well  for  shot  arrays  as  for  individual  explosions  Th e u ess 
of  this  discriminant  for  explosions  with  high  F values  suggests  that  the 
waves  from  explosions  are  not  due  to  tectonic  strain  release,  althoug  ^ 
conclusion  is  probably  not  well-established  without  better  aximut  a con 
The  preferred  discriminant  is  to  measure  the  maximum  S on  any  componen  , 
see  Table  IV  to  compute  and  discriminate  on  Ms  - MSHEAR. 

From  Figure  41a,  b we  can  also  see  that  the  ratio  SH/SV  also  is  not  a 
good  discriminant.  Figures  32  and  33  suggest  that  one  explanation  for  h 
is  a substantial  amount  of  site-t„-site  "rotation"  of  SV  into  SH;  possibly 
to  ray  refraction  out  of  the  great-circle  path. 

investigation  of  negative  discrimination  as  a technique  has  shown  that 
its  threshold  for  explosions,  assuming  about  0.01  probability  of  a false  Alarm, 
is  within  0.1-0.2  M magnitude  units  of  the  earthquake  threshold  • 

probability  of  detection  by  2 or  more  stations  for  a 25-station  world-wide 

network. 


,, 


Topics  for  future  research  include: 


• Investigation  of  long-period  P/LR  ratios; 

• Improved  distance-amplitude  relations  using  more  LPE  data; 

• Intensive  investigation  of  Asian  earthquakes  and  explosions 
using  NORSAR; 

• Investigation  of  such  puzzles  as:  Why  are  MILROW  and  CANNIKIN 

so  different.  Why  does  SJTX  have  such  a large  S wave  for 
CANNIKIN?  What  is  the  cause  of  the  large  SH  on  FI  and  E2  at 
LASA  for  CANNIKIN? 

• Better  azimuthal  coveragi  of  NTS  events  using  WWSSN  data. 
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A (degrees) 

Figure  2.  Travel  times  of  S to  thirty  degrees  as  read  at  North  American 
stations  for  the  Nevada  Test  Site  explosions  HALFBEAK  and  GREELEY.  The 
ordinates  scale  is  in  km.  From  Nuttli  (1969). 
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SH  AMPLITUDES 
GREELEY 


Figure  3.  SH  amplitudes  for  GREELEY  from  Nuttli  (1969).  The  smooth 
curve  through  the  data  has  been  drawn  by  hand.  The  shallow  sloping 
curve  through  the  peak  of  the  data  is  proportional  to  A-*?,  the  relation 
suggested  by  Evernden  (1969). 


diagrams  ol  mid-horizontal  section  through  plaster  cast 
by  small  explosion  in  soft  clay.  Illustrates  asymmetrvs 
From  Wright  and  Carpenter  (1962). 
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SM  AMPLITUDES  — GREELEY 


Figure  6a.  Polar  plot  of  the  SH  amplitudes  for  GREELEY . 
sense  of  motion  is  similar  to  that  for  HALFBEAK,  as  seen 
From  Nuttli  (1969) . 


SH  AMPLITUDES—  HALF-BTAK 


Polar  plot  of  the  SH  amplitudes  for  HALFBEAK.  The  dashed 
es  branch  AB  arrivals  at  the  shorter  distances  from  branch 
Note  the  reversal  of  the  sense  of  motion.  Nuttli  (1969) 
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Figure  8.  Subsurface  vertical  particle  velocity  record  for  GAS BUGGY 
From  Perret  (1970). 
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Figure  9.  Subsurface  vertical  displacement  records  for  GASBUGGY 
Perret  (1970). 
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Figure  10.  Subsurface  radial  particle  velocity  records  for  GASBUOr.Y 
From  Perret  (1972b). 
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Figure  15.  Schematic  elevation  of  SALMON  created  cavity  and  location 
of  instrument  station  for  STERLING.  From  Perret  (1968). 
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Figure  16. 
STERLING.  1 
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Figure  19.  Horizontal  radial  displacement  for  SALMON.  Pairs  of  traces 
show  agreement  between  single  and  double  integration  of  records  from  two 
different  instruments  and  illustrate  degree  of  data  reliability.  From 
Perret  (1967). 
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Figure  20.  Horizontal  tangential  displacement  from  SALMON.  Pairs  of 
traces  show  agreement  between  single  and  double  integration  of  records 
from  two  different  instruments  and  illustrate  degree  of  data  reliability 
From  Perret  (1967). 
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Figure  21.  Probability  of  detection  of  lou„-period  P at  LASA  using 
FKCOMB  for  events  in  the  Kurils  and  Japan  as  a function  of  m,  . 


LASA  "S"  IB  DEC  71  KAM 
START  TIME  FROM  LASASITE  F4 
352  22  10  11  1971 


Figure  23.  An  SH  signal  on  a beam  for  an  event  in  Kamchatka  which 
not  detected  by  FKCOMB, 
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Figure  25.  Distance-amplitude  •-’.lationship  for  long-period  S waves. 

Data  used  was  the  largest  component  of  all  long-period  S waves 
detected  at  LPE  stations  January-April  1972  from  events  of  depth  less 
than  60  km  and  with  LR  detected  at  three  or  more  LPE  stations.  Smooth 
curve  is  drawn  by  hand  using  mean  values  as  indicated.  Crosses  are  fron 
the  smooth  curve  through  Nuttli’s  (1969)  data  seen  in  Figure  3 normally 
0 , and  the  solid  dots  are  from  Gutenberg  and  Richters  (1956)  surfac 
relationship  seen  in  Figure  1. 
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Figure  27.  S-wave  detection  data  for  Kuril  Island-Kamchatka  area.  Fi 
Strauss  (1972).  Smooth  curve  is  fitted  cumulative  normal. 
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Ll  CANNIKAN  "P"  VERTICAL 
START  TIME  FROM  LASASITE  F4 
310  22  5 6 1971 
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Figure  29.  Vertical  beams  and  individual 
wave  at  LASA  from  CANNIKIN. 


traces  for  the  long-period  P 


e 33.  Transverse  beams  and  individual  traces  for  the  long-period  S 
at  LASA  from  CANNIKIN. 


Rn 

inn 

HI 

ill 

|fk|ti 

r 


'T*  * ' 


n™ 

mviii 

jA'JI 

■ ■'Ml  M 1 

f A TT 

j/j  , i\,  * 

& 

*wv  v yj*  y ^Vr 

U^/V  '"  s,J 

H 

In! 

irjyi 

Wl 

ISVll 

fflZwniiMR 

g§T|B 

IIYA  tti^w. 

p 

V? 

K1I  K 

III 

pmiuiw 

WnVinNRrV 

uu 

lu^iij 

Kvl 

iiWiHin 

Ml 

n msi 

LL  MILROW  9 OCT  69  S=222126  TRANSVERSE 
START  TIME  FROM  LASASITE  F4 
275  22  17  12  1969 
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Figure  35.  Transverse  beam  and  individual  traces  for  the  long-period 
S wave  at  LASA  from  MILROW. 
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PERCENT  PROBABILITY  Of  OETECTION  OF  LR  FOR  EARTHQUAKES 


Figure  42.  Probability  of  detection  of  S by  > 1,  2,  or  4 stations  out  of 
25  and  probability  of  a false  alarm  (no  detection,  resulting  in  a decision 
that  an  earthquake  is  an  explosion  by  means  of  negative  discrimination)  as 
a function  of  M for  a 25-station  worldwide  network  after  Romney  (1971)  for 
an  epicenter  in  Kamchatka.  Noise  and  signal  standard  deviation  is  0.3 
magnitude  units.  Arrows  indicate  points  referred  to  in  the  text. 
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Figure  43.  Threshold  magnitude  for  99%  probability  of  detection  of  S by 
1 or  more  stations  by  a worldwide  network  of  23  stations  with  mean  zero- 
noise  levels  of  19  mp  except  for  15  mp  for  LASA,  NORSAR,  and  ALFA,  a 
signal-to-noise  ratio  for  detection  of  r=1.5,  and  a standard  deviation  o 
signal  and  noise  of  0.3  magnitude  units.  The  star  indicates  the  epicent 
assumed  for  the  probability  of  detection  versus  M plots  such  as  Figure 


Figure  44,  Threshold  magnitude  for  90%  probability  of  detection  of  S by 
2 or  more  stations  by  a worldwide  network  of  25  stations  with  mean  zero-to 
noise  levels  of  19  mp  except  for  15  mp  for  LASA,  NORSAR,  and  ALPA,  a 
signal-to-noise  ratio  for  detection  of  r=1.5,  and  a standard  deviation  of 
signal  and  noise  of  0.3  magnitude  units.  The  star  indicates  the  epicenter 
assumed  for  the  probability  of  detection  versus  M plots  such  as  Figure  42 
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Figure  45.  Threshold  magnitude  for  90%  probability  of  detection  of  S by 
4 or  more  stations  by  a worldwide  network  oT  25  stations  with  mean  zero— to— peak 
noise  levels  of  19  mp  except  for  15  mp  for  LASA,  N'ORSAR,  and  AI.I’A,  a 
signal-to-noise  ratio  for  detection  of  r=1.5,  and  a standard  deviation  of 
signal  and  noise  of  0.3  magnitude  units.  The  star  indicates  the  epicenter 
assumed  for  the  probability  of  detection  versus  M plots  sucli  as  Figure  42. 


Figure  46.  Probability  of  detection  of  S by  > 
and  probability  of  a false  alarm  (no  detection 
that  an  earthquake  is  an  explosion  by  means  of 
a function  of  Mg  for  a 25-station  worldwide  ne 
an  epicenter  in  Kamchatka.  Noise  and  signal  si 


Figure  47.  Threshold  magnitude  for  99%  probability  of  detection  of  S by 
1 or  more  stations  by  a wor  dwide  network  of  25  station  with  mean 
zero-to-peak  levels  of  19  my  except  for  15  my  for  LAS A,  NORSAR,  and  ALPA 
a signal-to-noise  ratio  for  detection  of  r=1.5,  and  a standard  deviation 
of  signal  and  noise  of  0.1  magnitude  units.  The  star  indicates  the 
epicenter  assumed  for  the  probability  of  detection  versus  M plots  such 
as  Figure  42. 


igure  48.  Threshold  magnitude  for  90%  probability  of  detec 
2 or  more  stations  by  a worldwide  network  of  25  stations  witl 

Spa'  a“Ps!enalT  °V9  mM  eXCept  f°r  15  for  LASA 

ALFA,  a signal-to-noise  ratio  for  detection  of  r=1.5  and  a 

deviation  of  signal  and  noise  of  0.1  magni.ude  units.  The  si 
as  pS^S:  3SS  f°r  the  r,robability  of  detection  versU! 


Figure  49.  Threshold  magnitude  for  90%  probability  of  detection  of  S by 
4 or  more  stations  by  a worldwide  network  of  25  stations  with  mean 
zero-to-peak  noise  levels  of  19  mp  except  for  15  mp  for  LASA,  NORSAR,  and 
ALPA,  a signal-to-noise  atio  for  detection  of  r=1.5,  and  a standard 
deviation  of  signal  and  noise  of  0.1  magnitude  units.  The  star  indicates 
the  epicenter  assumed  for  the  probability  of  detection  versus  M plots  such 
as  Figure  42. 


PROB.  OF  DETECTION 


PROB  Of  DETECTION 


Figure  50.  Probability  of  detection  of  1.R  by  >_  1,  2,  or  4 stations,  out 
of  25  and  probability  of  a false  alarm  (no  detection  resulting  in  a decision 
that  an  earthquake  is  an  explosion  by  means  of  negative  discrimination)  as 
a function  of  Mg  for  a 25-station  worldwide  network  after  Romney  (1971)  for 
an  epicenter  in  Kamchatka.  Noise  and  signal  standard  deviation  is  0.5 
magnitude  units. 
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Figure  51.  Threshold  magnitude  for  99%  probability  of  detection  of  LR  by 
1 or  more  stations  by  a worldwide  network  of  25  stations  with  mean 
zero— to— peak  noise  levels  of  19  m)J  except  for  15  mp  for  LASA,  NORSAR,  and 
ALPA,  a signal-to-noise  ratio  for  detection  of  r=1.5,  and  a standard 
deviation  of  signal  and  noise  of  0.3  magnitude  units.  The  star  indicates 
the  epicenter  assumed  for  the  probability  of  detection  versus  plots  su 
as  Figure  42. 
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4 or  more  stations  by  a worldwide  network  of  25  stations  with  mean 
zero-to-peak  noise  levels  of  19  m|i  except  for  15  m|J  for  LASA,  NORSAR,  and 
ALPA,  a signal-to-noise  ratio  for  detection  of  r=1.5,  and  a standard 
deviation  of  signal  and  noise  of  0.3  magnitude  units.  The  star  indicates 

as  Figure  42. 


PERCENT  PROBABILITY 


Pf  * Po(0|  = 1-Pjjlc  1) 


Figure  54.  Probability  of  detection  of  LR  by  1,  2,  or  4 stations  out  of 
25  and  probability  of  a false  alarm  (no  detection  resulting  in  a decision 
that  an  earthquake  is  an  explosion  by  means  of  negative  discrimination) 


as  a function  of  M for  a 25-station  worldwide  network  after  Romney  (1971) 
for  an  epicenter  in  Kamchatka.  Noise  and  signal  standard  deviation  is  0.1 


magnitude  units. 
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Figure  55.  Threshold  magnitude  for  99%  probability  of  detection  of 
1 or  more  stations  by  a worldwide  network  of  25  stations  with  mean 
zero-to-peak  noise  levels  of  19  my  except  for  15  my  for  LASA,  NORSAR 
ALPA,  a signal-to-noise  ratio  for  detection  of  r=1.5,  and  a standard 
deviation  of  signal  and  noise  of  0.1  magnitude  units.  The  star  indii 
the  epicenter  assumed  for  the  probability  of  detection  versus  M plo 
as  Figure  42.  s 
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Figure  56.  Threshold  magnitude  for  90%  [ 

2 or  more  stations  by  a worldwide  network  of  25 
zero-to-peak  noise  levels  of  19  r; 

ALPA,  a signal-to-noise  ratio  for  detection 
deviation  of  signal  and  noise  of 
the  epicenter  assumed  for  the 
as  Figure  42. 
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Figure  57.  Threshold  magnitude  for  90%  probability  of  detection  nf  tr  k 
4 or  more  stations  by  a worldwide  network  of  detfcblon  of  LR  by 

as  Figure  42.  Probability  of  detection  versus  Ms  plots  such 
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